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CHAPTER  I 


INTRODUCTION 


All  fixed  wing  aircraft  generate  trailing  vortex  wakes  whose 
circulatory  strength  is  dependent  on  flight  conditions  and  aircraft 
geometry.  Large  heavy  airplanes  can  generate  extremely  strong  trailing 
vortex  systems  during  landing  approaches,  takeoffs,  or  other  low  speed 
flight  conditions.  Wetmore  and  Reeder^^  and  many  others  have  studied 
this  problem  as  it  applies  to  aircraft  terminal  operations,  and  it  is 
well  established  that  lighter  and  smaller  aircraft  can  be  severely 
affected  by  encounters  with  these  strong  vortex  wakes.  When  this  occurs 
the  aircraft  ordinarily  experiences  severe  rolling  or  pitching  moments, 
or  a  combination  of  both,  depending  on  the  aircraft  attitude  with 
respect  to  the  vortex  system.  These  moments  may  cause  loss  of  control 
or  structural  damage  to  the  aircraft 8  If  the  aircraft  is  centered 
between  the  trailing  vortices  of  a  larger  airplane  it  may  experience  a 
severe  down  draft.  Although  larger  aircraft  are  not  as  strongly 
affected  by  these  vortices,  large  rolling  moments  can  still  be  developed 
which  could  be  disasterous  during  low  altitude  flight. 

These  vortices,  being  in  a  real  fluid,  have  a  central  region  of 

(2) 

rotational  motion  defined  as  a  vortex  core.  McCormick  ,  Hoffman  and 
Joubert^,  Spreiter  and  Sacks and  others  have  studied  the  problem 
of  predicting  the  geometry  of  a  vortex.  As  a  result,  a  great  deal  of 
information  has  been  assembled  concerning  vortex  geometry  and  behavior. 

The  trailing  vortices  of  a  fixed-wing  aircraft  are  formed  when  the 

vortex  sheet  shed  from  the  wing  rolls  up  downstream  into  two  discreet 

(2) 

vortices  located  inboard  from  each  wing  tip.  As  noted  by  McCormick  , 


2 

and  many  others,  this  roll  up  is  quite  rapid  and  occurs  within  two  or 
three  chord  lengths  of  the  wing  trailing  edge,  The  two  vortices 
rotate  in  opposite  directions  but  with  the  same  circulatory  strength. 

For  a  wing  the  circulatory  strength  of  these  vortices  is  the  midspan 
value  of  the  bound  circulation  on  the  wing.  Figure  1  is  a  sketch  of 
the  roll-up  of  the  vortex  sheet  shed  from  a  wing  in  forward  flight. 

Each  vortex  induces  a  velocity  on  the  other  causing  the  vortex  system 
to  move  downward  at  a  constant  rate  until  the  system  decays  or  meets  a 
boundary  such  as  an  airport  runway,  at  which  time  the  two  vortices 
spread  out  laterally  along  the  surface  until  they  decay. 

A  rotor  also  generates  a  vortex  wake.  A  vortex  sheet  is  shed  from 

each  blade.  Each  sheet  rolls  up  rapidly  outboard  to  form  a  tip  vortex. 

The  inboard  sheet  never  coalesces  into  a  strong  vortex  system  and  from 

(5) 

Lehman’s  studies  of  rotor  flow  patterns,  it  appears  to  have  no  effect 
on  the  far  field  wake  configuration.  There  are  basically  three  types  of 
far  field  wake  configurations  which  the  tip  vortices  then  form,  all 
being  variations  of  a  helical  trail  whose  boundaries  are  defined  by  a 
skewed  elliptical  cylinder  Lehman’s  studies  of  rotor  flow  patterns 
clearly  show  the  characteristics  of  each  type.  Referring  to  Figure  2, 
the  first  pattern  is  characterized  by  nearly  vertical  helical  vortex 
trails,  which  occur  at  low  tip  speed  ratios.  The  second  type  of  pattern 
occurs  at  moderate  tip  speed  ratios  where  the  skew  angle  is  becoming 
large.  The  tip  vortices  combine  to  form  two  oppositely  rotating 
vortices  similar  to  a  fixed  wing  vortex  system.  The  final  configuration 
is  that  of  a  coarsely  pitched,  highly  skewed  helix.  It  generally  occurs 
at  high  tip  speed  ratios  and  large  skew  angles.  Figure  3  is  an 
illustration  of  the  similarities  between  the  fixed-wing  and  rotor  wakes. 
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The  purpose  of  this  study  is  to  estimate  the  rotor  wake  geometry 

and  velocity  field  downstream  of  the  rotor  for  the  particular  case 

where  the  rotor  wake  exhibits  strong  similarities  to  the  fixed-wing 

wake,  A  method  is  developed  by  which  the  helicopter  far  field  wake  can 

be  readily  compared  to  the  far  field  wake  generated  by  an  equivalent 

wing.  With  this  method  an  estimate  of  helicopter  wake  decay  can  be  made 

(2) 

based  on  work  done  by  McCormick,  Tangier,  and  Sherrieb  for  trailing 

vortices  of  conventional  aircraft « 

The  conventional  methods  of  predicting  rotor  wake  geometry  and 

velocity  fields  in  forward  flight  by  the  use  of  analytical  models  are 

reasonably  accurate  in  predicting  blade  loadings  and  time  varying  flow 

fields  in  the  vicinity  of  the  rotor*  Landgrebe^\  Heyson^^,  and 
ro\ 

Crimiv  ,  to  name  a  few,  have  done  significant  work  in  this  area* 

These  methods,  however,  are  not  applicable  far  downstream  from  the  rotor 
because  they  do  not  consider  real  fluid  effects  in  the  wake  which  must 
be  considered  in  predicting  the  geometry  of  the  wake,  They  also  involve 
time  consuming  computer  programs  whose  accuracy  is  unnecessary  for  far 
field  wake  predictions . 

An  early  study  using  a  smoke  visualization  technique  and  a  hot-wire 

anemometer  to  measure  and  locate  rotor  blade  tip  vortices  was  done  by 

(9) 

Simons,  Pacifico,  and  Jones  ,  Their  study  indicated  that  the  outboard 
shed  vortex  sheet  was  fully  rolled  up  within  one  rotor  radius.  This 
result  is  clearly  validated  by  pictorial  results  from  a  later  study  done 
by  Lehman  concerning  model  helicopter  rotors, 

Landgrebe^  and  Crimi^^  use  distorted  wake  models  to  calculate 
blade  loads  and  near  field  velocities.  Both  of  these  analytical  models 
have  shown  the  tendency  for  the  modeled  rotor  wake  to  begin  rolling  up 
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into  two  trailing  vortices  downstream  from  the  rotor.  Heyson's^^  model 
for  a  rotor  wake  is  an  undistorted  classical  model.  Although  the  wake 
geometry  is  fixed,  the  downstream  velocity  field  appears  to  be  similar 
to  a  fixed-wing  wake  velocity  field®  In  a  later  publication  with  most 
of  the  same  information,  Heyson^^  compares  measured  velocities  down¬ 
stream  from  the  rotor,  which  occur  between  the  center  of  the  two 
trailing  vortices  in  the  plane  of  the  vortices,  with  calculated 
velocities  based  on  a  model  for  a  fixed-wing  far  field  wake.  He 
obtained  reasonable  results,  but  the  model  presented  was  only  designed 
for  calculating  the  center  line  velocities,  Heyson's  measured  velocities 
behind  an  actual  large  scale  model  rotor  also  clearly  show  the 
similarities  between  the  rotor  and  fixed-wing  far  field  wake.  The  tip 

(9) 

vortex  systems  generated  by  Hey son  s  large  scale  model  always  rolled 
up  into  two  trailing  vortices  within  a  rotor  radii.  One  of  the 
conclusions  that  Heyson  made  from  the  results  of  his  study  was  that  the 
far  field  rotor  wake  could  be  modeled  more  accurately  by  a  wing  than  a 
rotor  if  the  rotor  lift  coefficient  is  sufficiently  high. 

There  have  been  some  studies  done  on  the  aircraft-helicopter  wake 

(12) 

interaction  problem;  one  such  study  was  done  by  Connor  and  O' Bryan v  . 
Their  study  consisted  basically  of  flight  tests  with  an  airplane  and 
helicopter.  The  results  of  their  study  indicated  that  helicopter  wakes 
presented  potential  hazards  to  aircraft  operating  in  the  vicinity  of  the 
helicopter.  They  suggested  that  a  separation  interval  of  about  one 
minute  between  helicopter  and  airplanes  would  be  sufficient  for  safe 
operation  of  the  airplane  for  their  particular  generation  of  helicopters. 
It  should  be  noted  however  that  their  report  is  over  a  decade  old. 

Today,  helicopters  with  high  payload  capacities  are  much  larger  than 
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they  were  ten  years  ago,  indicating  that  the  helicopter  vortex- airplane 

interaction  is  now  a  more  serious  problem, 

A  general  approach  to  the  problem  of  estimating  far  field  rotor 

wake  geometry  and  velocities  is  described  below.  Both  analytical 

methods  and  experimental  data  are  incorporated  in  this  study.  In 

(9) 

addition  to  HeysonTs  measured  data  for  a  large  scale  rotor,  measure¬ 
ments  were  made  behind  a  model  rotor  by  the  author. 

Several  different  analytical  approaches  will  follow.  First,  a 
gross  treatment  of  the  rotor  is  made,  whereby  the  rotor  disk  is 
considered  to  be  a  low  aspect  ratio  elliptic  wing,  Geometric 
considerations  are  then  presented  to  estimate  the  trailing  vortex 
strength.  Next,  a  fixed  wake  model  of  the  tip  vortices  is  numerically 
integrated  to  obtain  velocity  fields  for  various  operating  conditions. 
These  velocity  fields  are  also  used  to  estimate  the  effective  trailing 
strength.  Finally,  an  analysis  of  the  rotor  wake  circulatory  strength 
is  made  by  using  an  analytical  model  of  a  flapping  rotor  in  forward 
flight. 

With  the  results  from  the  above  analyses,  together  with  experi¬ 
mental  data,  the  relationship  between  the  rotor  and  an  equivalent  wing 
are  established.  This  relationship  can  then  be  used  in  conjunction 
with  Reference  2  to  model  the  helicopter  far  field  wake. 
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CHAPTER  II 
ANALYSIS 


A  Gross  Treatment  of  the  Rotor  Disk  as  a  Wing 

An  elliptic  wing  of  semi-span  R  with  a  flight  speed  of  V  would 
have  a  mid-span  bound  circulation  of  given  by 


r 

OG 


2L 

TTpVR 


(1) 


Assuming  a  thrusting  rotor  with  a  rotor  radius  of  R,  and  a  flight  speed 
of  V,  to  be  analogous  to  the  elliptic  wing,  the  total  circulatory 
strength  of  the  trailing  vortices  would  be: 


T  =  (2) 

00  TTpVR  v  ' 

The  circulation  at  the  core  radius  of  a  vortex  has  been  shown  by 

(2) 

McCormick  et  alo  to  be  linearly  proportional  to  the  total 
circulation  of  the  vortex* 

r(a)  -  k  (3) 

The  constant  of  proportionality  k,  is  dependent  on  the  nature  of  the 
vortex*  For  small  scale  wing  models  k  has  a  value  of  .716,  which  comes 
from  a  solution  of  the  Navier-Stokes  equations.  For  large  scale  wings 
and  real  aircraft  k  has  a  value  of  approximately  0.16*  This  value 
results  from  empirical  measurements  presented  in  Reference  2.  This 
reference  also  shows  that  the  circulation  of  large  scale  vortices 

(3) 

follows  a  logarithmic  distribution  as  presented  by  Hoffman  and  Joubert 
The  rotor  lift  coefficient  will  be  defined  as  the  rotor  thrust 


divided  by  the  free  stream  dynamic  pressure  and  disk  area. 
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L  l/2pV27TR2 


(4) 


(13) 

Analysis  of  Grow's  data  on  the  maximum  induced  tangential  velocity 
of  a  trailing  vortex  shows  that  the  maximum  velocity  is  proportional  to 
the  wing  lift  coefficient  and  free  stream  velocity.  The  constant  of 
proportionality  is  independent  of  aspect  ratio  for  aspect  ratios  greater 
than  4,  with  a  value  of  0.625*  For  smaller  aspect  ratios  the  constant 
increases.  At  an  aspect  ratio  of  2,  the  value  os  0.7.  The  rotor  disk 
has  a  circular  plan  form  with  an  aspect  ratio  of  4/tt.  No  data  is 
available  for  the  constant  of  proportionality  for  this  aspect  ratio ,  but 
it  is  assumed  that  the  constant  has  a  value  of  at  least  0.7.  Thus 
app  roximat ely 


v.  -7  CL  V  (5) 

Combining  equations  (4)  and  (5)  yields  a  relationship  for  the  maximum 
induced  velocity  of  the  rotor  trailing  vortex  system. 


v. 

i 


1.4 


CT  VT 


(6) 


Where  is  the  tip  speed,  and  is  the  rotor  thrust  coefficient 
defined  as; 


CT  = 


2 

pirR  VT 


]i  is  the  tip  speed  ratio  given  by: 


(7) 


with  equations  (2) , 


(8) 


(3),  and  (6),  an  expression  for  the  core  radius 


can  be  found. 
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a  =  r  (9) 

1  .  4Tr 

Thus  by  comparing  the  rotor  to  an  equivalent  fixed  wing,  an 
estimate  can  be  obtained  for  the  core  size  of  the  rotor’s  trailing 
vortices.  An  alternate  method  will  now  be  presented  for  accomplishing 
the  same  thing. 

Detailed  Examination  of  the  Rotor  Vortex  System 

A  rotor  in  forward  flight  generates  a  bound  circulation  on  the 
rotor  blades  which  varies  both  radially  and  azimuthally.  This  leads  to 
the  supposition  that  the  trailing  vortices  from  each  side  of  the  rotor 
are  not  of  equal  strength.  If  true,  a  more  accurate  measure  of  the 
circulatory  strength  is  necessary. 

A  rotor  sheds  a  trailing  vortex  from  each  rotor  blade  tip.  For 
moderate  tip  speed  ratios  these  tip  vortices  can  combine  downstream 
from  the  rotor  to  form  two  trailing  vortices  as  noted  earlier.  If  the 
strength  of  one  of  the  combined  trailing  vortices  is  denoted  as  T  9 

then  an  amplification  factor,  A„,  will  be  defined  as  the  ratio  of  V 

f  °° 

to  the  product  of  the  number  of  blades  and  the  strength  of  the  rolled- 
up  tip  vortex  produced  by  one  blade* 

r 

00 

A  =  —  (10) 

BrB 

Denoting  the  quantity  BT  as  T  ,  the  total  strength  for  B  blades, 

D  0 

equation  (10)  is  simply: 


r 

00 


(ii) 
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An  approximate  estimate  of  the  amplification  factor  can  be  made 
solely  on  the  basis  of  the  geometry  of  the  rotor  tip  path.  First,  it 
is  assumed  that  the  vortex  from  any  single  blade  is  coincident  with  the 
blade1 s  tip  path.  Now,  for  a  rotor  system  traveling  in  the  negative  x 
direction  with  velocity  V,  the  equations  describing  the  tip  path 
coordinates  are. 


x  =  R  cos  (cj)  +  oot)  ~  Vt 

(12) 

y  =  R  sin  (<p  +  oot) 

(13) 

where  <j>  is  the  angle  defining  the  rotor  blade  position  at  time  t  =  0, 
oo  is  the  angular  velocity  of  the  rotor,  and  R  is  the  rotor  radius. 

Figure  4  is  a  plot  of  the  tip  path  of  one  rotor  blade.  If  the  coordinate 
system  origin  corresponds  to  time  t  -  0,  then  negative  time  denotes  a 
positive  x  coordinate^  The  substitution  \p  =  ~  oot  is  made  in  equations 
(12)  and  (13).  Non-dimensionalizing  by  the  rotor  radius  yields, 

x/R  =  cos  (( p  -  ip)  +  iJJUp  (14) 

y/R  =  sin  (<p  -  ip)  (15) 

For  any  x  station  there  can  be  from  one  to  an  infinite  number  of  ip 
values  which  satisfy  equation  (14),  depending  on  the  value  of  the  other 
parameters.  Each  solution  represents  a  separate  vortex  filament 
located  at  that  particular  x  station0  If  each  trailing  vortex  is 
comprised  of  one  half  of  the  filaments  located  at  any  x  station  which 
combine  to  form  that  vortex,  then  an  estimate  of  the  amplification 
factor  would  be  one  half  of  the  number  of  solutions  for  ip  in  equation 
(14).  Because  equation  (14)  is  a  transcendental  equation  a  general 
solution  is  difficult  to  obtain*,  An  approximate  solution  can  be  made, 
however.  We  begin  by  rearranging  equation  (14)  as, 


x/R  -  yijj  =  cos  ( (p  -  ip) 


(16) 
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Plotting  both  sides  of  equation  (16)  separately,  and  superimposing  them 
on  the  same  plot,  results  in  Figure  5.  It  is  evident  that  all  of  the 
solutions  to  equation  (16)  will  occur  within  the  region  defined  by 
|x/R  -  y^|  <;  10  The  length  of  this  region  is  then 

AijJ  =  2/y  (17) 

If  the  tip  speed  ratio  is  sufficiently  small,  the  number  of  inter¬ 
sections  of  the  cosine  function  with  the  ip  axis  and  the  number  of 
solutions  to  equation  (16)  will  be  approximately  equal*  The  amplifi¬ 
cation  factor  is  half  the  number  of  solutions,  therefore, 

Af  -  1/741  (18) 

The  important  behavior  to  note  from  this  simple  analysis  is  that  the 
amplification  factor  varies  inversely  with  tip  speed  ratio*  A  more 
accurate  estimate  of  the  constant  of  proportionality  will  be  made  later. 
This  approach  also  fails  to  consider  the  possibility  of  unequal 
strengths  for  the  two  trailing  vortices.  A  more  sophisticated  approach 
is  necessary 0 

Figure  6  is  a  sketch  of  the  tip  vortices  from  a  2-bladed  rotor  for 
three  separate  tip  speed  ratios  *  As  the  tip  speed  ratio  increases  it 
is  evident  that  the  advancing  side  has  a  greater  number  of  filaments 

(14) 

that  can  combine  to  form  a  trailing  vortex.  Eisenhuth  and  McCormick 
noted  that  when  two  vortex  filaments  intersected  each  other  there  was  a 
tendency  for  the  filaments  to  combine  into  separate  loops  0 

Figure  7  is  an  illustration  of  this  behavior  for  the  rotor  vortex 
geometry  of  Figure  6*  This  figure  further  illustrates  the  uneven 
distribution  of  filaments  as  noted  above.  As  the  tip  speed  ratio 
increases  the  advancing  side  gains  a  greater  portion  of  the  shed  tip 
vortices „  This  indicates  that  the  advancing  side  of  the  rotor  would 
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Figure  7.  Tip  Vortices  of  Figure  6  Combined  into  Separate  Filaments. 
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have  the  stronger  trailing  vortex  strength  if  the  tip  vortices  are  of 
constant  strength.  Thus,  as  the  tip  speed  ratio  increases  equation  (18) 
may  become  progressively  inaccurate,  since  it  is  based  on  assuming  that 
the  two  trailing  vortices  are  of  equal  strength. 


Flapping  Rotor  in  Forward  Flight 

A  more  accurate  estimate  for  the  amplification  factor  can  be  made 
using  the  relationships  derived  for  a  flapping  rotor  without  taper  or 
twist  by  Gessow  and  Myers^~^0  The  lift  per  blade  integrated  over  one 
period  of  revolution  is: 

L*?  p*oVEVl2  [1+f  ^  +  If1  (19) 

0 

where  A  is  defined  as  the  inflow  ratio,  0q  is  the  hub  pitch  angle,  and 

a  is  the  slope  of  the  lift  curve  of  the  rotor  airfoil  section.  The 
o 

bound  circulation  on  any  rotor  blade  is  found  to  be: 

^  3. 

r(x,iM  =  ca  e  VT[x  +  0“  +  \  v*  0^  +  sintKy-  0^-  x) 

0  o  o 


-  cos  ip  (y  +  X  -gb  ]  (20) 

0  o 

where  Bq,  a^,  b^,  are  respectively  the  coning  angle,  longitudinal 
flapping  and  lateral  flapping.  Integrating  over  one  period  of 
revolution  to  find  the  mean  bound  circulation  at  any  station  on  the 
blade  yields: 


F(x) 


1 

2tt 


2tt 


r(x,i|0dTp  =  fcao0oVT 


(21) 


0 

From  a  moment  balance  on  the  rotor  blade  the  longitudinal  flapping 


coefficient  is: 
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1-1 /2u2 


<t°o 


+  2  A) 


(22) 


Assuming  that  the  maximum  bound  circulation  occurs  at  the  tip  of  each 
rotor  blade,  then  each  shed  tip  vortex  has  a  mean  circulation  of 


rB-KVl!1+f +2  I11 

0  0 


(23) 


Combining  equations  (19)  and  (23)  yields: 


L  =  ieRVTrB 


,  ,  3  2  ,  3  A 

i  +  2  y  +  2  e 


1  +  9~+  2  y  91 
0  0 


(24) 


for  the  resultant  lift  per  cycle  per  blade  *  For  B  blades  the  rotor 
lift  then  becomes: 


L=?EV° : :  r~r  i 

1+0  +  2  y  0 
o  o 


(25) 


The  lift  for  a  wing  of  semi -span  R  and  flight  speed  V  with  an  elliptical 
circulation  distribution  is, 


L  =  f  pVRT^  (26) 

where  is  the  mid-span  circulation  of  the  wing.  The  ratio  of  V V  to 
rQ>  i.e.  the  amplification  factor,  is  then, 


A 


f 


2 

3ttu 


♦! 


y  + 


A_ 

0 


(27) 


0  0 

Equation  (27)  is  plotted  in  Figure  8  versus  tip  speed  ratio  for  various 
A/G^'s.  For  a  model  described  in  Reference  16,  A/ 0q  was  found  to  be 
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approximately  -.152  for  a  tip  speed  ratio  of  .31.  From  Figure  8  it  can 


be  seen  that  for  this  case  the  amplification  factor  behaves  essentially 

2 


as : 


Af  3  iTy 


(28) 


with  some  slight  change  in  magnitude.  It  should  be  noted  here  that  the 
variables  involved  in  equation  (27)  are  interdependent  and  must  therefore 
be  determined  for  a  given  flight  condition  to  obtain  an  accurate  esti¬ 
mation  of  the  amplification  factor.  For  a  first  estimate  the  term  in 
the  brackets  of  equation  (27)  can  probably  be  assumed  equal  to  unity. 

Equation  (27)  is  still  based  on  the  assumption  that  the  two 
trailing  vortices  for  the  rotor  system  have  equal  strengths.  This 
assumption  is  reasonable  because  a  trimmed  rotor  must  behave  in  the 
same  manner  as  a  fixed  wing;  the  induced  downwash  should  be  symmetrical 
about  the  ^  =  0,  180°  line.  This  occurs  when  the  trailing  vortices 
have  equal  strengths.  But,  Figure  7  indicates  that  this  is  not  true  if 
one  assumes  a  mean  circulation  strength  for  the  tip  vortices  as  is  done 
to  obtain  equation  (27),  The  following  analysis  will  show  why  the 
trailing  vortices  do  have  approximately  equal  strengths. 

The  tip  circulation  from  equation  (20)  is: 


rB(i|>)  =  J  cao  e0vT  I1!  +  r2  cos1^  +  r3  sinl(J] 


(29) 


where  and  r  ,  can  be  found  from  equation  (20).  It  is  assumed 

that  this  is  also  the  maximum  bound  circulation  and  therefore  the  tip 
vortex  circulation.  From  a  moment  balance  of  the  rotor  in  Reference  15, 
the  lateral  flapping  coefficient  b^  is 


4yB 


b.  = 


1  2 

3(1+  ) 


(30) 
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From  equation  (22)  and (30)  F^,  and  become 


p  ri.5  2.  A  ....  .  1  2.  ,  /ft  1  2. 

^  =  [1  +  T  |J  +  q“  (1  +  J  U  )]/(l  "  J  P  ) 

°0 

3  2 

6  7  +  ~  y 

2  3  0  \  .  1  2; 

o  1  +  2^ 


(31) 


(32) 


„  ,5,1  2  2A.  1  ,2. 

r3  y  b  +  2  y  0  ^  ^  2  y  ^ 

o 


(33) 


Equation  (29)  can  be  rewritten  as: 


r* 


rB<« 


Kwi 


r  r 

2  3 

1  +  —  cos  \\)  4-  -p—  sin  \p 

1  1 


(34) 


From  equations  (31),  (32),  (33),  and  (34),  it  is  apparent  that  as  the 
tip  speed  ratio  increases  the  advancing  sideTs  circulation  decreases 
while  the  retreating  sidefs  circulation  increases.  On  the  other  hand. 
Figure  7  shows  that  the  number  of  vortices  on  the  advancing  side 
increases  with  y  while  the  number  on  the  retreating  side  decreases. 

Thus  Figure  7  leads  to  the  conclusion  that  the  two  trailing  vortices  of 
the  rotor  system  could  have  nearly  equal  strengths  and  that  the 
amplification  factor  as  defined  by  equation  (27)  will  be  valid. 

To  evaluate  the  constants  on  the  sine  and  cosine  functions  of 
equation  (34),  the  inflow  ratio  and  coning  angle  must  be  known.  The 
coning  angle  can  be  found  from  the  following  equation  obtained  by 
Gessow  and  Myers from  trim  considerations. 


*0  = 


1  A 

-rcpa  0^ 

2  o  o 


/ 3 A  ,  1  , .  ,  2. 
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(35) 


23 


where  I  is  the  blade  moment  of  inertia,  M  the  weight  moment  of  the 
1  w 

blade,  and  0)  is  the  angular  velocity  of  the  rotor  blade.  The  inflow 
ratio  can  be  found  from  the  following  expression: 


a 


T 


1 

2  2  2 

y(u  +  xV 


(36) 


where  a  is  the  angle  of  attack  of  the  rotor  disk  plane*  The  thrust 
coefficient  is  simply, 


„  1  q  rl  1  2  L  1  A  - 

CT  2  aoa0o  [3  +  2  y  +  2  0  1 

0 


(37) 


where  O  is  the  rotor  solidity. 


a  = 


BL 

ttR 


(38) 


The  longitudinal  flapping  coefficient  can  also  be  found  from  helicopter 
trim  considerations^0^ 


“l  +  “  -  'f 

Equations  (22),  (36),  (37),  and  (39),  define  a  system  of  equations  of 
five  unknowns;  C^,  a,  a^,  A,  0^,  Given  the  helicopter  geometry,  tip 
speed,  flight  speed,  body  weight,  and  drag  characteristics,  and  a 
desired  thrust  coefficient,  the  above  set  of  equations  can  be  solved  for 
the  unknown  quantities*  Once  these  are  known,  the  coning  angle  and 
lateral  flapping  coefficient  can  also  be  found* 

Figure  9  is  a  plot  of  equation  (34)  for  a  fictitious  helicopter 
operating  at  two  different  tip  speed  ratios.  The  helicopter 


characteristics  are  assumed  to  be: 
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W  =  7000  lbs. 

VT  =  650  fps 
p  =  . 00237  slug/ft3 
0  =  .  06 


a  =  5.73  C  /rad 
o 

2 

f  =  16  ft  (equivalent  flat  plate  area) 
CT  =  ,00555 
R  =  20  ft. 
y  -  o  1;  .2 


M 


Y 


•  00228 

4 

cpa  R 
o 


5.5 


This  figure  clearly  illustrates  that  the  tip  vortices  shed  from  each 
blade  are  stronger  on  the  retreating  side  than  those  on  the  advancing 
side  and  that  the  difference  becomes  greater  as  the  tip  speed  ratio 
increases. 


Calculation  of  the  Far  Field  Velocity  Using  a  Potential  Model 

To  see  how  accurately  equation  (28)  models  the  amplification 
factor,  a  simple  potential  model  of  a  helicopter  wake  is  used  to 
calculate  the  velocity  field.  This  calculated  velocity  field  can  then 
be  compared  to  a  fixed  wing  velocity  field  to  calculate  the  theoretical 
amplification  factor. 

For  simplicity  in  calculating  the  rotor  far  field  velocities  it  is 
assumed  that  the  tip  vortex  from  each  blade  is  coincident  with  the  tip 
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path  for  all  time,  rather  than  using  a  free  wake  analysis*  The  major 
difference  between  the  two  models  is  the  relative  location  of  the  rotor 
wakes  in  the  vertical  direction*  This  fixed  wake  method  allows  a 
considerable  savings  in  computer  time  for  calculations  of  the  far  field 
velocity  of  the  rotor  vortex  systems. 

For  the  fixed  wake  model  the  Biot-Savart  Law  can  be  used  to  find 


the  velocity  at  any  point  in  the  x-y  plane »  The  relationship  is 

p  p  x  ds 


dv 


4ttR  ~  i  3 
P 


(40) 


where  p  and  s  are  non-dimensional  vectors  based  on  the  rotor  radius. 

For  this  model  p  is  the  vector  from  the  point  of  interest  to  any  point 
of  the  vortex  filament,  and  ds  is  the  differential  length  of  the  vortex 
filament  at  that  point.  Therefore, 


p  «  (x  -  X  )  I  +  (y  -  y  )  j  (41) 

and 

ds  =  d  x  I  +  dy  j  (42) 


where  x  and  y  are  defined  by  equations  (14)  and  (15)  respectively  and 
have  been  non-dimensionalized  by  the  rotor  radius.  Using  the  Biot-Savart 
Law  for  a  system  where  the  limits  on  ip  are  infinite,  the  velocity  at  any 
point  in  the  x-y  plane  becomes: 


VvV  =  m 


f°°  r (4>)  { (x-x  )  cos  (<M0+(y-y  )  (sin(<M0+U)  1 


{(x-x  )2  +  (y-y  )2} 
o  o 


3/2 


dip  (43) 


For  simplicity  we  assume  that  the  circulation  strength  differs  only 
slightly  from  the  mean  value.  From  equation  (34)  it  was  shown  that  this 
assumption  increases  in  error  as  the  tip  speed  ratio  increases,  but  for 
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a  first  approximation  it  is  assumed  to  be  sufficiently  accurate.  Using 
a  constant  circulation  in  equation  (43)  results  in  the  following 
expression  for  the  non-dimensional  induced  velocity. 


4ttR  /  v  r* 
p  V. (x  ,y  )  =  -k 

1  1  O  O 


f°  {(x-x  )cos(t()-^)+(y-y  )  (sin(<|)-i(j)+u)  } 


r  /  s2,.  .2,3/2 

{(x-xq)  +(y-yQ)  } 


dip 

(44) 


The  limits  of  integration  have  been  taken  over  all  ijj  because  this  model 
is  used  to  obtain  the  amplification  factor  in  the  ultimate  wake.  This 
equation  was  evaluated  by  numerical  integration. 

A  wing  in  forward  flight  generates  a  vortex  sheet  which  rolls  up 
into  two  vortices  which  lie  inboard  from  the  wing  tips.  The  far  field 
velocities  can  be  calculated  by  assuming  that  the  tip  vortices  are  two 
discreet  vortex  filaments  starting  at  the  wing  tips  and  continuing  to 
infinity.  However  the  inboard  movement  of  the  trailing  vortices  will 
be  ignored  in  this  analysis  since  it  is  assumed  that  the  rotor-generated 
far  field  vortices  also  have  no  inboard  movement. 

The  Biot-Savart  Law  for  a  straight  line  filament  of  finite  length 
reduces  to 

v(h)  =  ^~(cos  y  +  cos  3)  (45) 

where  y,  3,  and  h  are  defined  in  Figure  10.  For  a  wing  semi-span  of  R 
centered  in  the  x-y  plane  origin,  the  far  field  velocity  can  be  found 
by  the  use  of  equation  (45)  applied  to  the  two  filament  system.  For 
this  trailing  vortex  system  defined  in  the  x-y  plane,  the  non-dimensional 
induced  velocity  in  the  z  direction  becomes: 


4itR  ,  , 

~r~  v(V 


where  y^  has  been  non-dimensionalized  by  the  semi -span. 


(46) 
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As  stated  previously  the  analytical  model  presented  here  is  based 
on  a  fixed-wake  geometry.  Because  of  this  the  velocity  field  at  any  x 
station  is  highly  dependent  on  the  blade  position  at  time  t  =  0,  To 
overcome  this  dependency  it  is  desirable  to  take  results  from  various  x 
stations  defined  over  one  region  of  periodicity,  For  each  x  station  an 
amplification  factor  can  be  found  when  the  rotor  velocity  field  is 
compared  to  the  equivalent  fixed  wing  velocity  field.  This  is  done  by 
calculating  the  necessary  circulation  of  the  wing  such  that  the  velocity 
fields  described  by  equations  (44)  and  (46)  match,  From  these  results 
amplification  factors  for  each  x  station  within  one  region  of  periodicity 
are  calculated  and  then  averaged  to  obtain  a  time  averaged  value. 

The  region  of  periodicity  is  found  in  the  following  manner.  This 
region  is  the  change  in  the  downstream  coordinate,  x,  as  the  rotor  blade 
goes  through  one  revolution.  The  period  of  one  revolution  is 


With  the  use  of  equation  (14)  the  change  in  x  is  found  to  be: 

Ax  -  2TTRy  (48) 

This  region  can  be  broken  into  several  small  intervals  as  described 
above  to  evaluate  the  amplification  factor. 

Numerical  integration  of  equation  (44)  was  done  on  an  IBM  360/67 
digital  computer.  Since  the  limits  of  integration  extend  over  all  tjj  a 
finite  integration  angle  must  be  used,  A  description  of  the  truncation 
of  ijj  is  given  in  Appendix  II,  The  resulting  velocity  fields  were 
compared  with  the  velocity  described  by  equation  (46)  for  a  fixed  wing, 
and  an  amplification  factor  was  found.  Figure  11  presents  the  predicted 
amplification  factor  as  a  function  of  tip  speed  ratio  for  a  single 
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bladed  rotor  and  a  six  bladed  rotor.  The  amplification  factor  was 
found  to  be  essentially  independent  of  the  number  of  rotor  blades  used 
and  given  by: 


A 


f 


25 

U 


(49) 


This  agrees  reasonably  well  with  equation  (28)  and  appears  to  justify 
the  gross  treatment  of  the  rotor  as  an  equivalent  wing. 
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CHAPTER  III 
EXPERIMENTAL  STUDY 


Description  of  Experimental  Apparatus 

The  model  used  for  the  velocity  field  measurements  was  a  two-bladed 
untwisted  teetering  rotor  with  a  radius  of  ll!l  and  a  hub  radius  of  1.5". 
The  rotor  section  was  a  NACA  0012  airfoil  section  with  a  chord  of  2"  • 

The  rotor  hub  was  made  of  steel  and  the  rotor  blades  of  aluminum.  The 
model  rotor  was  dynamically  balanced  for  rotational  speeds  of  up  to 
1800  rpm,  the  maximum  speed  that  the  motor  could  turn  the  rotor.  All 
measurements  were  taken  with  a  hub  pitch  angle  of  10  degrees.  The 
rotor  was  driven  by  a  variable  speed  3/4  hp  motor.  The  rotor  was 
tested  in  the  4x5  foot  subsonic  wind  tunnel  at  The  Pennsylvania  State 
University,  Department  of  Aerospace  Engineering.  For  all  measurements 
the  rotor  disk  plane  was  aligned  with  the  free  stream  velocity  such 
that  the  rotor  disk  had  no  angle  of  attack* 

Three  different  tip  speed  ratios  were  tested.  For  each  case  the 
velocity  measurements  downstream  from  the  rotor  model  were  taken  with  a 
45°  hot-wire  probe  mounted  on  a  motor-driven  traversing  mechanism  and 
used  in  conjunction  with  a  battery  operated  anemometer.  The  signal  is 
converted  to  a  linearly  proportional  signal  using  King’s  Law.  (See 
Appendix  I) .  Figure  12  is  a  diagram  of  the  analog  circuitry  used  for 
the  conversion*  A  EAI  model  TR-20  Analog  Computer  was  used  for  this 
purpose.  The  signal  is  first  filtered  to  remove  high  frequency 
information  and  then  passed  into  the  analog  computer.  The  resulting 
signal  is  then  sent  into  a  EAI  model  1130  Variplotter.  The  y-plot 
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signal  is  driven  by  a  voltage  from  a  potentiometer  mounted  on  the 
traversing  mechanism  and  the  velocity  signal  is  recorded  on  the  Vari- 
plotter.  Figure  13  is  a  block  diagram  of  the  electronic  equipment  used 
in  this  set-up,  and  Figure  14  is  a  sketch  of  the  experimental  set-up. 

Test  of  the  Hot-Wire  Anemometer 

In  order  to  test  the  validity  of  the  hot-wire  as  a  measuring 
device  and  to  gain  confidence  in  the  hot-wire  output  as  compared  to 
known  velocity  fields  measured  with  a  vortex  meter,  the  hot-wire  was 
first  used  to  measure  the  velocity  field  behind  a  fixed  wing  for  three 
different  operating  conditions 0 

The  wing  used  had  a  NACA  0015  section  profile,  a  semi-span  of  18", 

and  a  chord  of  6" 0  Figure  15  is  a  plot  of  the  maximum  measured  induced 

tangential  velocity  versus  the  wing  C  „  Good  agreement  with  the 

JL 

(2) 

results  of  McCormick  et  ale  was  obtained  for  the  maximum  induced 
velocity o  Figure  16  is  a  typical  plot  of  the  induced  velocity  field 
versus  radial  position  for  the  measurements  taken,,  The  induced  velo¬ 
city  has  been  non-dimens ionalized  by  the  free  stream  velocity „  The 

core  radius  for  this  case  was  found  to  be  approximately  0o47  inches . 

2 

This  gives  a  circulation  at  the  core  radius  of  3.57  ft  /sec.  The 

2 

estimated  mid-span  circulation  for  this  core  was  5o43  ft  /sec*  The 
value  for  k  in  equation  3  for  this  case  was  found  to  be  CL  658,  All 
measurements  were  taken  33"  downstream  of  the  rotor  trailing  edgeQ 
From  these  results  it  appears  that  the  trailing  vortices  behave  in  the 
manner  predicted  by  the  Navier-S tokes  solution  as  expressed  in  Reference 
2.  The  results  also  indicated  that  the  hot-wire  would  be  a  reliable 
measuring  device. 


Schematic  Diagram  of  Experimental  Apparatus  in  Test  Section. 


NACA  0015 


b  =  18" 
c  =  6" 

V  =  29  fps 
x  =  39"  from  L.E 


I - 1 — 

-4.0  -3.0 


Figure  16. 
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Measurements  of  Rotor  Velocity  Field  and  Sample  Calculations 

Having  gained  confidence  from  the  fixed  wing  tests,  the  hot-wire 
was  then  used  to  measure  the  velocity  field  behind  the  teetering  rotor 
model  as  described  previously*  Three  different  operating  conditions 
were  measured.  In  each  case  the  velocity  field  measured  along  the 
y-axis  downstream  of  the  rotor  was  recorded  on  the  Variplotter «  This 
was  done  at  intervals  of  1  inch  in  the  z  direction*  All  velocities 
were  measured  at  x/R  =  3.55*  In  this  manner  the  vortex  cores  were 
located  for  both  advancing  and  retreating  sides 0  Because  the  single 
45°  hot-wire  probe  is  only  calibrated  for  velocities  in  the  plane 
defined  by  the  intersection  of  the  mean  tunnel  velocity  and  the  wire 
itself,  only  calibrated  signals  are  obtained  when  the  probe  passes 
through  or  close  to  the  vortex  core  center.  These  core  centered  measure¬ 
ments  are  the  measurements  used  in  this  study* 

As  noted  earlier,  the  rotor  disk  plane  was  aligned  with  the  free 
stream  velocity.  Lateral  flapping  was  negligible  for  all  cases  studied. 
For  a  zero  angle  of  attack  of  the  rotor  disk  plane,  the  inflow  ratio  is 


where  w  is  the  downwash  of  the  rotor  a  Assuming  a  constant  downwash 
equal  to  the  mean  downwash  of  the  rotor  in  forward  flight,  and  taking 
finite  hub  radius  and  tip  losses  into  consideration,  equation  (37) 
becomes 
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where  T  has  the  value  .97  and  h  is  the  hub  radius  non-dimens ionali zed 

D 

by  the  rotor  radius. 

The  far  field  circulation  can  then  be  found  by  the  use  of  equation 
(2)  rewritten  in  terms  of  the  thrust  coefficient  and  tip  speed  ratio. 


r 

00 


2  CL 


The  bound  circulation  at  the  tip  is  found  from  Equation  (20).  For  this 
teetering  rotor  model  it  is: 


rw  =  i  caoe0vT[1 


2y0, 


+  y  sin  i(j] 


Table  1  is  a  listing  of  the  operating  parameters  used  for  the 
three  different  cases  studied,  and  the  calculated  parameters  of 
interest*  Sample  calculations  for  one  case  are  presented  in  Appendix 
III. 
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CHAPTER  IV 
RESULTS 


Figure  17  is  a  typical  example  of  the  velocity  profile  measured  by 
the  hot-wire  anemometer  behind  the  rotor  as  it  was  plotted  by  the 
Variplotter  for  this  study*  Figure  18  is  the  same  plot  with  the  small 
scale  oscillations  faired  out,  and  with  velocity  magnitude  added*  These 
measurements  were  taken  at  the  z  station  representing  the  vortex  center 
for  which  the  core  diameter  and  maximum  tangential  velocities  were 
largest . 

Because  the  measurements  taken  for  this  study  were  for  a  small 
scale  model,  Heyson’s^"^  measured  velocities  behind  a  large  scale 
model  were  incorporated  in  this  study  to  obtain  more  information 
concerning  rotor  far  field  velocities 0  Heyson’s  measurements  were 
taken  at  various  x,  y,  and  z  stations  underneath  and  behind  a  rotor 
model.  The  rotor  used  was  a  teetering  type  with  untapered,  untwisted 
blades.  The  rotor  blades  had  a  NACA  0012  airfoil  section*  The  rotor 
radius  was  7,5  feet  with  a  solidity  of  .0534o  In  order  that  Heyson’s 
data  correspond  as  accurately  as  possible  to  data  taken  in  this  study, 
the  velocity  measurements  at  the  x/R  station  of  3*14  in  Heyson’s  study 
were  used,  these  measurements  being  the  closest  to  the  x/R  station  of 
3.55  used  in  this  project’s  study*  All  velocity  measurements  for 
Heyson’s  study  were  made  with  pitot  tube  rakes*  Because  of  this,  the 
velocity  measurements  are  time  averaged  values*  These  measurements  were 
made  in  the  NASA  Langley’s  full  scale  wind  tunnel.  For  each  case  the 
measurements  at  the  z  station  which  had  the  greatest  tangential 


V 


Advancing  Side 
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velocities  at  the  core  were  usech  The  information  obtained  from 
Hey son's  study  is  listed  in  Table  20 

Figure  19  is  a  plot  of  core  radius  versus  tip  speed  ratio  for  both 
the  data  taken  by  Heyson^^  and  that  obtained  by  the  author.  Equation 
(9)  predicts  that  the  core  radius  should  be  independent  of  tip  speed 
ratio,  and  only  a  function  of  whether  the  vortex  system  is  laminar  or 
turbulent.  The  results  of  Figure  19  suggest  that  this  is  the  case 
since  the  larger  rotor  exhibits  the  smaller  relative  core  size.  However, 
there  is  a  tendency  for  the  core  radius  to  decrease  in  size  for  both 
rotors  as  the  tip  speed  ratio  increases 0 

Figure  20  is  a  plot  of  the  constant  k  defined  by  equation  (3)  for 
the  experimental  data  obtained  based  on  the  circulation  measured  at 
the  core  radius  and  the  estimated  total  circulation  of  the  rotor 0  As 
can  be  seen  from  Figure  20  the  smaller  model  resulted  in  higher  k  values 
than  those  obtained  from  Heyson's  measurements*  In  both  cases  the 
derived  k  values  appear  to  be  high,  particularly  in  view  of  the  fact 
that  the  values  for  the  small  model  are  greater  than  unity* 

j! 

Figure  21  is  a  plot  of  the  measured  maximum  induced  tangential 
velocities  divided  by  the  thrust  coefficient  and  tip  speed  versus  tip 
speed  ratio  Equation  (6)  is  also  plotted  on  this  figure.  The 
relationship  as  expressed  by  equation  (6)  is. 


1,4 

y 


The  least  squares  fit  of  the  data  to  the  parameters  involved  in 
equation  (6)  yielded  the  following  relationship. 
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Table  1.  Data  Taken  at  The  Pennsylvania  State  University. 


Case  No. 
and  Side 

y 

VT 

(fps) 

V 

(fps) 

V.* 

i 

(fps) 

a 

(in) 

CT 

r(a) 

(f2ps) 

r 

00 

(f2ps) 

1  Adv. 

.1285 

144.0 

18.5 

16.9 

3.6 

.0112 

31.86 

23.06 

1  Ret. 

15.5 

3.2 

25.97 

2  Adv. 

.201 

129.5 

26.0 

12.3 

3.6 

.0134 

23.19 

15.81 

2  Ret. 

13,5 

3.2 

22.62 

3  Adv. 

.301 

86.4 

26.0 

11.5 

3.0 

.0157 

L8.06 

8.285 

3  Adv. 

10.2 

0.4 

2.14 

*  Averaged  Over  Inboard  and  Outboard  Side  of  Each  Core 


Table  2.  Data  Taken  by  Heyson,  Reference  11. 


Case  No. 
and  Side 

y 

VT 

(fps) 

V 

(fps) 

V.  * 

X 

(fps) 

a 

(ft) 

CT 

T(a) 

(f2ps) 

✓“N 

03 

8  CMa 

1  Adv. 

.095 

500 

47.5 

34.1 

.855 

.0032 

183.0 

253.0 

1  Ret. 

23.0 

1.434 

207.0 

2  Adv. 

.140 

500 

o 

o 

27.4 

.645 

.0037 

110.5 

199.0 

2  Ret. 

15.7 

1.575 

155.4 

3  Adv. 

.232 

450 

104.5 

12.7 

.645 

.0032 

51.4 

93.3 

3  Ret. 

13.4 

1.15 

96.3 

*  Averaged  Over  Inboard  and  Outboard  Side  of  Each  Core 


This  is  close  to  the  prediction  of  equation  (6),  even  though  the 
equation  was  based  on  a  wing  with  an  aspect  ratio  of  2. 
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CHAPTER  VI 

CONCLUSIONS  AND  DISCUSSIONS 


From  the  analyses  and  empirical  measurements  performed  in  this 
study,  it  appears  that  the  far  field  wake  of  a  rotor  can  be  adequately 
described  by  treating  the  rotor  as  a  circular  wing.  As  a  result,  the 
following  relationships  were  found  to  be  valid: 

1.  The  maximum  induced  tangential  velocity  about  a  trailing 
vortex  from  a  rotor  is  proportional  to  the  product  of  the  thrust  coef¬ 
ficient  and  tip  speed  divided  by  the  tip  speed  ratio.  The  constant  of 
proportionality  is  approximately  1.48,  The  maximum  tangential  velocity 
is  essentially  the  same  as  that  which  one  would  obtain  by  treating  the 
rotor  as  a  low  aspect  ratio  fixed  wing, 

2.  The  size  of  the  core  radii  of  the  trailing  vortices  of  the 
rotor  system  is  proportional  to  the  rotor  radius,  and  appears  to  be 
independent  of  other  operating  conditions «  The  constant  of  proportion¬ 
ality  is  a  function  of  the  nature  of  the  vortex  system,  i.e.  whether 
the  flow  in  the  vortex  is  laminar  or  turbulent.  This  behavior  for  core 
radii  size  is  similar  to  the  fixed  wing  core  radius  size  behavior  which 
is  proportional  to  the  mid-span  chord  of  the  wing,  in  that  the  rotor 
radius  is  proportional  to  the  mean  chord  of  a  circular  wing. 

3.  The  large  and  small  scale  rotors  exhibit  scale  effects  with 
regard  to  the  circulation  at  the  core  radius.  It  appears  that  the 
trailing  vortex  geometry  for  small  scale  rotors  can  be  modeled  by  decay 
laws  based  on  the  Navier-S tokes  equations  while  the  trailing  vortex  of 
the  large  scale  rotor  is  more  accurately  modeled  by  a  logarithmic 


circulation  model  incorporation  empirical  decay  laws  for  fixed 
wing  aircraft  as  described  by  reference  2. 


54 


REFERENCES 


1.  Wetmore,  J.  W,  ,  and  Reeder,  J0  P.,  "Aircraft  Vortex  Wakes  in 

Relation  to  Terminal  Operations,"  NASA  TND-1777,  April,  1963. 

2.  McCormick,  B.  W. ,  Tangier,  J.  Le ,  and  Sherrieb,  Ha  E, ,  "Structure 

of  Trailing  Vortices,"  Journal  of  Aircraft,  Vol.  5,  No.  3,  May- 
June,  19680 

3.  Hoffmann,  E.  R, ,  and  Joubert,  P0  N0,  "Turbulent  Line  Vortices," 

Journal  of  Fluid  Mechanics,  Vole  16,  pp0  395-411,  1963c 

4.  Spreiter,  J*  R0 ,  and  Sacks,  A0  He,  "The  Rolling  Up  of  the  Trailing 

Vortex  Sheet  and  its  Effect  on  the  Downwash  Behind  Wings,"  Journal 
of  Aeronautical  Sciences,  Vol0  18,  No0  1,  ppa  21-23,  January,  1951. 

5.  Lehman,  A.  F. ,  "Model  Studies  of  Helicopter  Rotor  Flow  Patterns," 

USAAVLABS  Technical  Report  68-17,  April,  1968. 

6.  Landgrebe,  A.  J.,  "An  Analytical  Method  for  Predicting  Rotor  Wake 

Geometry,"  AIAA  Paper  #  69-196. 

7.  Heyson,  H.  H@ ,  "Analysis  and  Comparison  with  Theory  of  the  Flow 

Field  Measurements  Near  a  Lifting  Rotor  in  the  Langley  Full  Scale 
Tunnel,"  NACA  TN  3691 * 

8.  Crimi,  P.,  "Prediction  of  Rotor  Wake  Flows,"  CAL/USAAVLABS 

Symposium  Proceedings,  Vol.  I  Propeller  and  Rotor  Aerodynamics, 
June,  1966. 

9.  Simons,  1.  A0 ,  Pacifico,  R0  Eo ,  and  Jones,  J.  P.,  "The  Movement, 

Structure  and  Breakdown  of  Trailing  Vortices  from  a  Rotor  Blade," 
CAL/USAAVLABS  Symposium  Proceedings,  Vol.  I  Propeller  and  Rotor 
Aerodynamics,  June,  196 6 „ 

10.  Lehman,  A*  J. ,  "Test  Section  Size  and  Influence  on  Model  Helicopter 

Rotor  Performance,"  USAAVLABS  Technical  Report  71-6,  March,  1971. 

11.  Heyson,  H.  H. ,  "Induced  Velocities  Near  a  Lifting  Rotor  with 

Nonuniform  Disk  Loading,"  NACA  TR  1319. 

12.  Connor,  A.  B.,  and  0TBryan,  "A  Brief  Evaluation  of  Helicopter  Wake 

as  a  Potential  Operational  Hazard  to  Aircraft,"  NASA  TND-1227, 
March,  1962 o 

13.  Grow,  T.  Lo ,  "The  Effect  of  Wing  Geometry  and  Lower  Surface  Boundary 

Layer  on  the  Rolled-Up  Tip  Vortex,"  M0So  Thesis,  The  Pennsylvania 
State  University,  1967 . 


55 


14.  Eisenhuth,  J.  J.,  and  McCormick,  B.  W. ,  "An  Observation  on  the 

Vortex  System  of  Dual-Rotation  Propellers,"  Journal  of  Aeronautical 
Sciences,  Vol.  20,  No.  6,  September,  1953. 

15.  Gessow,  Alfred,  and  Myers,  G.  C. ,  Jr.,  Aerodynamics  of  the  Helicopter. 

New  York:  Fredrick  Ungar  Publishing  Co.,  1952. 

16.  McCormick,  B.  W.  Aerodynamics  of  V/Stol  Flight.  New  York:  Academic 

Press,  1967. 

17.  PISA  Electronics  and  Service  Manual  for  Type  55005  Battery  Operated 

CTA.  DISA  Electronics  A/S,  Herlev,  Denmark. 


56 


APPENDIX  I 


Basic  Principle  of  the  Constant  Temperature  Hot-Wire  Anemometer 
The  hot-wire  anemometer  works  on  the  principle  that  a  thin 
electrically  heated  wire,  suspended  between  two  prongs,  conducts  heat 
to  the  fluid  surroundings®  The  heat  transfer  from  the  wire  obeys  the 
two-dimensional  heat  transfer  relationship  for  a  cylinder  as  described 
by  L.  V.  King®  The  formula,  known  as  King's  Law,  states  that  the 
rate  of  heat  loss  of  the  wire,  Q,  is  related  to  the  fluid  velocity  U 
by  the  following  relationship: 

1/2 

Q  =*  (T  -  Tq)  (A1  +  B1  u  ) 

where  T  is  the  wire  operating  temperature,  T^  is  the  fluid  temperature, 
and  A^,  are  empirical  constants. 

The  rate  of  heat  loss  must  equal  the  power  generated,  P,  in  the 
wire  by  the  electric  current,  I,  when  the  wire  is  held  in  thermal 
equilibrium.  That  is, 

2 

P  =  I  R;  T  -  T  =  constant 
o 

where  R  is  the  wire  operating  resistance® 

Under  thermal  equilibrium  conditions 

R  =  R  (1  +  a  (T  -  T  )  ) 
o  o 
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where  is  the  resistance  of  the  wire  at  fluid  temperature  and  a  is  a 
constant  of  proportionality.  Thus,  King's  Law  can  be  rewritten  as 


~-A2  +  B2 


R-R 


This  relationship  has  been  found  to  be  valid  in  atmospheric  air  at 
velocities  in  the  range  of  .5  m/sec  to  100  m/secs. 

A  feedback  system  keeps  the  wire  temperature  constant  with  the 
output  signal  being  a  voltage  E  which  obey's  King's  Law  in  the  following 
form: 


E 


2 


A  +  B 


u1/2 


The  hot-wire  is  sensitive  only  to  velocities  normal  to  the  wire* 

If  the  wire  axis  is  inclined  from  the  direction  of  the  flow  velocity  by 
an  angle  0,  then  the  relationship  becomes: 


where 


e2  = 


A  +  B  u 


,1/2 


u'  =  u  sin  0 


For  a  two-dimensional  flow  u'  becomes 


u'  =  V  sin  0  +  u  cos  0 

if  V  is  a  constant*  then  u  can  be  easily  measured  with  the  hot-wire. 
The  final  relationship  for  u  as  a  function  of  E  and  the  empirical 
constants  A  and  B  is  then: 
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U  =  CSC  0 


f[^]2  - 


V  sin  0} 


Figure  22  illustrates  the  geometry  of  the  hot-wire  probe  as  applied  to 
a  two-dimensional  flow  with  one  vector  component  constant. 


Figure  22.  Hot-Wire  and  Two-Dimensional  Flow  Geometry 
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APPENDIX  II 


Truncation  of  ^ 

Because  it  is  necessary  to  integrate  equation  (38)  numerically  with 
a  digital  computer,  a  method  for  defining  the  length  of  vortex  system 
necessary  to  obtain  solutions  approaching  solutions  of  the  infinite 
system  is  required.  The  following  method  of  defining  a  satisfactory, 
finite  length  was  used. 

Consider  an  infinite  vortex  filament  with  constant  strength.  From 
the  Biot-Savart  Law  the  velocity  at  any  radial  distance  from  the 
filament  is 


v(h) 


r  2 

4ttR  h 


where  h  =  — 


For  a  finite  length  filament  the  velocity  is  just: 


F  cos  y  +  cos  6 
4ttR  h 


The  error  in  the  velocity  calculation  due  to  using  a  finite  length  in 
place  of  an  infinite  length  filament  would  be 

cos  y  +  cos  6 
E  -  1  -  2 

For  simplicity  consider  only  the  case  where  the  limits  on  the 
finite  length  filament  are  symmetrical  about  the  point  of  interest, 
then  cos  y  =  cos  8,  therefore 


E  =  1  -  cos  y 
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From  Figure  23  it  can  be  seen  that: 


cos  y  = 


£ 


-7 


-  h 


2 


where  £  is  a  non-dimensional  length.  This  yields: 


E  = 


1  - 


£ 


Solving  for  £  as  a  function  of  h  and  the  desired  error  yields: 


£  =  h 1L 
4 


2  2 
2E  -  E 


For  the  symmetrical  case  the  total  segment  length  must  be  L  =  2£.  For 
a  desired  accuracy  this  segment  length  is  then  used  for  the  limiting 
length  on  the  numerical  integration.  This  length,  £,  in  turn  can  be 
directly  related  to  the  angle  ifj  with  the  use  of  equation  (14)  to 
obtain  an  integration  limit  on  \p0  For  this  study  an  error  of  5  %  was 
deemed  acceptable. 


Figure  23.  Geometry  fi 

Cu 
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APPENDIX  III 


Sample  Data  Reduction 

From  Table  1,  the  tip  speed  ratio  is 

y  =  V/VT  =  18.5/144.0  =  .1285 

For  the  NACA  0012  airfoil  the  slope  of  the  lift  curve  is  approximately 

a  -  5.73  C„ /rad 
o  Jt 

The  solidity  for  this  rotor  is 

C  =  Bc/ttR  =  (2)  (2") /tt(11")  =  .11575 
The  non-dimensional  hub  radius  is 

h  =  1. 5"/ll. 0"  =  .1364 

The  hub  collective  pitch  angle  was  10°  or  .174  radians. 

With  this  information,  the  thrust  coefficient  for  the  teetering 
rotor  can  be  found  approximately  by  the  following  expression: 


Cm  = 


ir  aa 
2  o  0 


3~(tb2  -  h3> +  i  <tb  - h) 


1  + 


2  2 
tb  -  h 
8 


a  a 
o 


or  for  this  particular  case. 


01122 
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Therefore,  the  total  circulation  for  the  far  field  vortes  should  be 

ra  =  2  CT  RVT/y  =  23.06  ft2/sec. 

On  the  advancing  side  the  following  measured  values  were  obtained: 

v,  =  16*9  fps 
i 

a  =  3.6" 

Therefore,  the  circulation  at  the  core  radius  is 

T(a)  ■  2'rrav^  =  31.86  ft2/sec. 

On  the  retreating  side  the  following  measured  values  were  obtained: 

v.  =  15.5  fps 
i 

a  =  3.2" 

Therefore,  the  circulation  at  the  core  radius  is 

r(a)  =  25.97  ft2/ sec. 

The  k  values  for  this  case  are  1.38  and  1.13  for  the  advancing  and 
retreating  sides  respectively. 


